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ABSTRACT 

Context. Massive protostars have associated bipolar outflows with velocities of hundreds of km s _1 . Such outflows can 
produce strong shocks when interact with the ambient medium leading to regions of non-thermal radio emission. 
Aims. We aim at exploring under which conditions relativistic particles are accelerated at the terminal shocks of the 
protostellar jets and can produce significant gamma-ray emission. 

Methods. We estimate the conditions necessary for particle acceleration up to very high energies and gamma-ray 
production in the non-thermal hot spots of jets associated with massive protostars embedded in dense molecular 
clouds. 

Results. We show that relativistic Bremsstrahlung and proton-proton collisions can make molecular clouds with massive 
young stellar objects detectable by the Fermi satellite at MeV-GeV energies and by Cherenkov telescope arrays in the 
GeV-TeV range. 

Conclusions. Gamma-ray astronomy can be used to probe the physical conditions in star forming regions and particle 
acceleration processes in the complex environment of massive molecular clouds. 
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1. Introduction 

Massive stars are formed in the dense cores of massive cold 
clouds (Garay & Lizano 1999, and references therein). The 
accumulation of gas in the core might proceed through pre- 
vious stages of fragmentation and coalescence with the pro- 
gressive result of a massive protostar that then accretes 
from the environment (e.g. Bonnell et al. 1997, Stahler et 
al. 2000) or through direct accretion onto a central object 
' of very high mass (e.g. Rodriguez et al. 2008 -RMF08-; 
see Shu et al. 1987 for the basic mechanism). In any case, 
the prestellar core is expected to have angular momentum, 
which would lead to the formation of an accretion disk. 
The strong magnetic fields inside the cloud that thread the 
disk should be pulled toward the protostar and twisted by 
the rotation giving rise to a magnetic tower, with the con- 
sequent outflows, as shown by numerical simulations (e.g. 
Banerjee & Pudritz 2006, 2007). 

Evidence of molecular outflows is found through 
methanol masers, which are likely associated with shocks 
formed by the interaction with the external medium (e.g. 
Plambeck & Menten 1990). However, the most important 
evidence for outflows comes from the detection of thermal 
radio jets. These jets are observed to propagate through 
the cloud material along distances of a fraction of a par- 
sec (e.g. Marti et al. 1993 -MRR93-). At the end point 
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of the jets, hot spots due to the terminal shocks are ob- 
served in several sources. In a few cases, these hot spots are 
clearly non-thermal, indicating the presence of relativistic 
electrons that produce synchrotron radiation (e.g. Araudo 
et al. 2007 -ARA07-, 2008). 

A population of relativistic electrons in the complex en- 
vironment of the massive molecular cloud in which the pro- 
tostar is being formed will produce high-energy radiation 
through a variety of processes: inverse Compton (IC) scat- 
tering of infrared (IR) photons from the cloud, relativis- 
tic Bremsstrahlung, and, if protons are accelerated at the 
shock as well, inelastic proton-proton (pp) collisions. If such 
radiation is detectable, gamma-ray astronomy can be used 
to shed light on the star forming process, the protostar 
environment, and cosmic ray acceleration inside molecular 
clouds. 

This work is devoted to discuss under what conditions 
the terminal shocks of jets from massive protostars can ef- 
ficiently accelerate particles, and produce gamma rays that 
may be detectable by the Fermi satellite and Cherenkov 
telescopes in the near future. The model developed for 
the calculations is essentially different from the phenome- 
logic model presented by ARA07, since now the dynamics 
of the jet termination region is characterized, the shock 
power estimated, the conditions for particle acceleration 
analyzed, and the role of hydrodynamical instabilities for 
non-thermal radiation explored. In short, the acceleration 
and emission are consistently modeled together with the 
hydrodynamics in a more physical scenario. Our study is 



2 



V. Bosch-Ramon et al.: Massive protostars as gamma-ray sources 



partially based on early works on non-thermal emission in 
young stellar objects (YSO), as those by Crusius-Watzel 
(1990) and Hcnriksen et al. (1991), but we develop further 
some important aspects of the hydrodynamics-radiation re- 
lation, and focus on massive YSOs and the feasibility of 
their detection with the present observational facilities. 



2. Physical scenario 

A massive YSO, or a group of them, are deeply embed- 
ded into a molecular cloud. The protostar heats the cloud 
in such a way that it can be detected as a strong IR 
source, with luminosities in the range Lir ~ 10 4-5 L 
~ if) 38-39 erg s" 1 , whereas the optical counterpart is ob- 
scured by the cloud. Masses and sizes of these clouds are of 
the order of ~ 10 3 M Q and few pc, respectively (e.g. Garay 
& Lizano 1999), and the densities in the regions in which the 
massive YSOs are found typically span n c ~ 10 3 — 10 6 cm" 1 
(see Araudo et al. 2008 and references therein). 

As already mentioned, collimated outflows are present 
in most of massive YSOs, and thermal radiation have been 
detected up to distances of 10 16 — 10 18 cm from the central 
star. These jets have temperatures of ~ 10 4 K and move at 
speeds (vj) from several hundreds to ~ 1000 km s" 1 (e.g. 
MRR93, Marti, Rodriguez and Reipurth 1995 -MRR95-). 
This shows that they are strongly supersonic flows with 
Mach numbers M <; 10. The kinetic luminosities of these 
jets are expected to be L } ~ 10 36 erg s _1 (e.g. MRR95, 
ARA07). 

In some cases (see Araudo et al. 2008), non-thermal ra- 
dio lobes have been detected at distances of Z- } ~ pc, with 
sizes of - 1% of this distance (MRR93, Garay et al. 2003 
-GAR03-). Given the directions, sizes and distances from 
the core, the lobes are compatible with being produced by 
the head of the jet. These radio lobes are likely generated by 
the strong terminal shocks of the jets, which also ionize the 
shocked material. Magnetic fields should also be present, 
since they play an important role supporting the cloud be- 
fore the gravitational collapse, allowing the required high 
densities in the cores to be achieved (e.g. McKee & Ostriker 
2007) . Estimates from Zeeman measurements of the cloud 
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magnetic field give values B c ~ 0.3n c g mG (e.g. Crutcher 
1999), where n C 5 = n c /(10 5 cm -3 ) is the cloud density. 
Under these conditions, particles can be accelerated up to 
relativistic energies via diffusive shock (Fermi I) acceler- 
ation (DSA, e.g. Bell 1978; see Drury 1983 for a review). 
These particles would produce the radiation of non-thermal 
nature found in the lobes, and could generate significant 
emission in a broad spectral range, from radio to gamma 
rays. 

Some amount of thermal ultraviolet (UV)/X-ray pho- 
tons is expected from the shocked material. This radiation 
will suffer strong photo-electric absorption in the regions of 
the cloud next to the jet head and will ionize the surround- 
ing material improving the conditions for efficient particle 
acceleration (e.g. Drury et al. 1996). On the other hand, 
fast radiative cooling of the shocked material can affect the 
lobe dynamics, and also reduce the efficiency of DSA, but 
could increase the detectability of massive YSO at high 
energies because of the density enhancement. Finally, ion- 
ization losses of radio emitting electrons and free-free ab- 
sorption by the ionized medium could affect significantly 




Fig. 1. Sketch of the termination region of the jet of a 
massive YSO. Two shocks of different strengths and veloc- 
ities, depending on the jet-medium properties, will form. 
Electrons and protons can be accelerated in the shocks, 
and generate non-thermal emission via interaction with the 
ambient matter, magnetic and radiation fields. The shocked 
material will also produce thermal radiation. 

the radio spectrum. In some sources, free- free emission may 
dominate the radio band. 

In Fig. [1] we sketch the scenario in which non-thermal 
emission is produced through the acceleration of electrons 
and protons in the jet termination regions. 

3. On the physical nature of the lobes 

We assume that the non-thermal radio lobes are the regions 
in which the protostellar jets terminate. The action of the 
jet head on the external medium leads to two shocks, one 
moving in the cloud material and another one in the jet 
itself; these are the bow shock and the reverse shock, re- 
spectively. These shocks would be the accelerators of the 
relativistic particles generating the observed non-thermal 
radio emission. 

3.1. Dynamics of the jet termination shocks 

An important parameter determining the shock character- 
istics is the jet (j) to cloud (c) density ratio x = n- 3 /n c . For 
fiducial values of the jet properties, say v- } ~ 10 8 cm s _1 , 
Lj ~ 10 36 erg s" 1 and jet radius R s ~ 10 16 - 10 17 cm (as- 
suming that lobe and jet radii are similar), we obtain jet 
densities in the range rij ~ 10 2 — 10 4 cm" 3 at the location 
of the lobe; then, \ ~ 10~ 4 - 10. 

The value of Xi together with Vj, determines the speed 
of the bow shock (e.g. Blondin et al. 1989): 

v hs « (1 + X~ 17 TS » (0-01 - 0.8)«j 

w 10 6 - 10 s v js cm s" 1 , (1) 

where v- 3 s = Wj/(10 8 cm s" 1 ). The reverse shock velocity 
is v r ~ Vj — 3/4 i>b s . The life time of the jet, which can be 
written as tuf e « Zj / Vbs > can be also expressed as a function 
of and the parameters Z pc — Zj/(1 pc), Xo.i = x/0-l ; 

tUfa » Z- s /v hs - 3 x 10 10 Z pc v" 1 s , when x > 1 , (2) 

and 

tlife ~ Zj/v bs rs 10 11 Z pc Vjg 1 Xai /2 s > wnen X < 1 • (3) 
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Equations [Hand [3] show that the jet advance takes place 
in two different regimes depending on the jet-medium den- 
sity ratio, which depends in turn on the source age. As 
long as the jet lateral pressure is larger than that of the 
surrounding medium, the jet expands freely and thereby 
Uj oc R? oc Z?, i.e. the jet has a conical shape. For very 
young jets, when % > 1, the advancing jet head speed is 
constant (t>b s ~ v j): an d therefore the jet lifetime is simply 
t oc Z- } . However, the jet head gets diluted. Since x ^j~ 2 i 
at some stage \ < 1 and t becomes oc Z?. This Zy-t de- 
pendence implies that most of the sources will be observed 
when x < 1- 

At some point, the jet expansion is stopped by the exter- 
nal pressure and the jet density becomes roughly constant. 
When it happens depends on the shocked jet material pres- 
sure away from the reverse shock, but it is expected that 
X <C 1. If values as low as x ~ 10~ 4 are reached the bow 
shock will move with a speed below the Alfven speed (i>a) 
in the cloud and will not be a shock anymore. 

The luminosity crossing the shock surface in the jet and 
the cloud (reverse and bow shock, respectively) is: 



(7r/2)i??m p n jiC i; r 3 bs 



3 x 10 35 R 



16 '*J,c 



3 «r,bs8 er S S 



(4) 



where 



■*j,c3 



= 7i jiC /(10 3 cm 3 ) and v^ hsa 



« r .bs/(10 8 cm s" 1 ). The sum of L s from both shocks 
cannot be larger than Lj . 

When x < 1, the reverse shock to bow shock luminosity 
ratio is ~ X -1 ^ 2 - This fact, together with the Zj-t depen- 
dence mentioned above, implies that most of the time the 
reverse shock will be more powerful than the bow shock. On 
the other hand, for x 3> 1 the reverse shock will be very 
weak and the bow shock radiation very faint. Therefore, 
lobe detection is favored when the jet has expanded enough 
to reach x Si 1, being the reverse shock the best place for 
particle acceleration unless x ~ 1, when both shocks have 
similar properties. 

3.2. The postshock regions 

For purely adiabatic jet-medium interactions, the formation 
of the reverse shock is accompanied by a strong widening of 
the jet termination region as in extragalactic FRII sources 
(e.g. Kaiser & Alexander 1997) and possibly also in micro- 
quasars (e.g. Bordas et al. 2009). However, the conditions 
in massive YSOs are likely different. As noted by Blondin 
et al. (1989), if the cooling length Z coo i of the shocked ma- 
terial in either the reverse, the bow shock, or in both, is 
< .Rj, the jet head will not expand significantly. Using the 
cooling function C(T) - 7 x icr 19 !? 1 - - 6 (e.g. Myasnikov 
et al. 1998) at the temperatures relevant here, adopting 
a density and speed downstream right after the shock of 
4nj jC and t | r ,bs/4, respectively (strong, non-relativistic and 
lowly magnetized shock), and assuming full ionization in 
that region, ^ coo i is: 



cool 



3kTv Tths 
32n hc C(T) 



10 



18 ,,4.2 
,bs8 



n j,c3 Cm - 



(5) 



Despite a specific-case treatment and detailed simulations 
would be required to characterize the fine evolution of the 
jet head, there is a wide range of realistic parameters for 



which the bow shock, and even the reverse shock, can be 
radiative. This is compatible with the relatively small size 
of non-thermal radio lobes (see Sect. [2]). In case one or 
both shocks are not radiative, the material cools through 
adiabatic expansion farther than Rj from the shock. The 
adiabatic cooling timescale is similar to the advection (or 
escape) time of the material in the downstream region: 



to 



4 i?j/w r ,bs = 4 x 10 s R il6 v r l s8 s . 



(6) 



Because of radiative cooling, the compression ratio A, 
or downstream to upstream density ratio, will grow with 
the distance downstream from A = 4 to ~ 3 x 10 3, u 2 8 , or 

to ~ 20 v S 8 n X J C 3 BZ_\ ■ The former upper limit for A comes 
from the fact that the temperature stops falling around 
T ~ 10 4 K, and the latter from the enhancement of the 
magnetic field (B = 10~ 3 £?_3 G) pressure, which limits 
the compression (see Blondin et al. 1989). 

When x < 1) the material downstream is 1/x to 
Abs/xA r times denser in the bow shock than in the re- 
verse shock. Since a force is exerted by the downstream 
material of the reverse shock on that of the bow shock, 
Ray leigh- Taylor (RT) instabilities in the contact disconti- 
nuity between both shocks can develop. This phenomenon 
will distort the jet head on a timescale: 



= t RT ~ Rj/x 



'1/2, 



o -i nQ 7-1 I— 111 —1 

3 X 10 %6Xo.001 U j8 s ; 



(7) 



where Xo.ooi = (X-^-r/Ab s )/0.001 is the effective density 
jump in the contact discontinuity. Therefore, for i^fc > 
£rt, complex shock structures should form (Blondin et al. 
1989) accompanied by strong mixing of material from both 
shocks. This mixing can lead to an increase of the effective 
density downstream the reverse shock. To account for this, 
we have simply parameterized the density downstream the 
reverse shock as 4nj < AFrij < Ab s n- C , where F is a free 
parameter that says how much the density departs from the 
adiabatic value, in y Downstream the bow shock, mixing 
should not affect significantly the density, and hence F ~ 1. 

Concerning the dynamic role of the jet and cloud mag- 
netic fields at the jet termination region, Bj and B c , re- 
spectively, we will assume hereafter that they are negligible. 
This applies as long as B hC <C B eq , where: 



£ Oq = 30L^ 2 ^ j8 1/2 mG 



(8) 



is the magnetic field of equipartition between the magnetic 
and the jet kinetic energy density. If B hC B cq , it is also 
expected that the magnetic field in the shock regions, B, 
should be well below -BccJB 

4. Particle acceleration in the lobes 

4.1. Acceleration and cooling processes 

The energies and luminosities that non-thermal particles 
may achieve depend on the efficiency of particle acceler- 
ation, which depends in turn on the lobe properties. As 
noted above, DSA can operate in the fast strong shocks 
at the termination of the jets accelerating particles up to 
relativistic energies. For a parallel, lowly magnetized (i.e. 



Note that B cq relates to the jet kinetic energy density, not 
to the one of the (unshocked) cloud. In the latter, the magnetic 
field may be indeed dynamically more important. 
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v r ,bs 3> va upstream), non-relativistic strong shock, in the 
test particle approximation and with diffusion coefficient 
D, the acceleration rate is (e.g. Protheroe 1999): 

fp 3 ,7-1 f V ^s\ 2 o 

^e,pgain — « V~ J 9 

« 1.5 x 10~ 5 cT 1 v 2 bs8 B_ 3 GeV s" 1 , (9) 

where e and p stand for electrons and protons, respectively, 
and d = D/Db, with Db = cr g /3 being the diffusion coef- 
ficient in the Bohm limit (r g = E/qB is the gyroradius of 
a particle with energy E). The acceleration timescale can 
be written as: 

tgain = E/E e . p gain « 6.7 X f 4 d« r ; b 2 sg SZ^GeV S , (10) 

where E GcV = J5/(lGeV). 

At the maximum energy of particles, t ga i n becomes equal 
to the shortest timescales among synchrotron, IC (Thomson 
regime) and relativistic Brcmsstrahlung losses for electrons 
(see Blumenthal & Gould 1970), inelastic pp collision losses 
for protons (see Kelner et al. 2006), and diffusive particle 
escape and jet life time for both electrons and protons. The 
relevant loss timescales, t\ oss = —E/E\ oas , are given by the 



following expressions: 

t sync ra 4 x 10 11 BZj E G l v s , (11) 

tic « 1-6 x 10 13 u^_ 9 E G l Y s , (12) 

tBrems « 3.5 X 10 10 U^ 3 FZq 1 S , (13) 

t pp ^5xW w n^ s F 1 - 1 s, (14) 

t diff » 1.5 x 10 12 tT 1 R 2 16 B_3 £ G cV s , (15) 

tufa « 10 n Zp C Ubis, (16) 



where uir-9 = uir/10 9 erg cm 3 is the energy density of 
the IR photon field and F 10 = F/10 (recall iq =0.1 in the 
bow shock). Concerning the densities, the compression of 
the material due to the shock is already accounted. Under 
the adopted assumptions, particles spend on average ~ 4 
times longer in the downstream than in the upstream re- 
gion, and the conditions downstream the shock will deter- 
mine the particle evolution. Adiabatic cooling is not con- 
sidered in the accelerating region. 

From the timescales presented above, we obtain the fol- 
lowing maximum energies depending on the dominant en- 



ergy loss mechanism and source age: 

E m ax sync « 2.4 x 10 3 cT 1/2 « r , bs8 BZ\ /2 GeV , (17) 

E max ic « 1.5 x 10 4 d- 1 ' 2 v r , hs8 uVlil B 1 / 2 GeV , (18) 

£max Brcms « 5.2 X 10 5 dT 1 < bs8 B_ 3 ^3 F^ 1 GeV , (19) 

E max P p « 7.5 x 10 5 cT 1 « 2 bs8 S_ 3 nr^3 i^ 1 GeV , (20) 

E max diff w 4.7 x 10 3 cT 1 u r , bs8 B_ 3 i? jl6 GeV , (21) 

E max i if e w 1.5 x 10 5 gT 1 Z pc Wl , bs8 B_ 3 GeV , (22) 



This shows that, if d ~ 1, electrons and protons can reach 
very high energies. It turns out that synchrotron and diffu- 
sive escape are the dominant mechanisms to limit accelera- 
tion, although relativistic Bremsstrahlung and pp collisions 
can dominate for large densities. We notice that in radia- 
tive shocks, Rj should be substituted by Z coo i in Eqs. ([16]) 
and (122]) . 

We note that the accelerated proton to electron num- 
ber ratio a cannot be established from first principles. It is 
considered here as a phenomenologic parameter. 

4.2. Required conditions for efficient DSA 

In DSA particles are scattered by magnetic inhomogeneities 
back to the shock before escaping downstream/upstream. 
These inhomogeneities should not move faster than the 
shock itself, since otherwise the particles would not be 
affected by the shock velocity jump. If it were the case, 
stochastic Fermi II particle acceleration could be at work 
(Fermi 1949). However, in that case most of the jet kinetic 
energy will not be available and little energy will go to non- 
thermal particles. Therefore, efficient particle acceleration 
requires the shock to be super-Alfvenic in the upstream 
region, i.e. va <C w r ,bs, where: 



t, A « ^4^-^ « 2.1x10^-^ cm s 4 , (23) 

and Xi = 0.1 Vjo.i is the ratio of ionized to total density. 
Thus, X; should not be too low. 

Downstream the shock, the material is likely ionized 
through particle collisions, and the magnetic field is ex- 
pected to be quite disordered. Upstream the shock, UV/X- 
ray radiation produced in the shock region can photo- 
ionize the medium unless «bs 10 7 cm s _1 (or \ 
0.01 fjg 2 ), since then photons cannot ionize hydrogen. From 
the recombination and photo-ionization timescales, t Te ~ 
10 11 XZq\ nZ c3 s and t p h ~ 10 7 R 2 6 v 2 s s, respectively 
(the recombination rate and ionization cross section are 
given in Seaton 1959 and Morrison & McCammon 1983; 
-L35 corresponds here to the ionizing photon field luminos- 
ity), it seems likely that upstream the shock t pn < t rc , and 
therefore X; — ► 1. In the radiative regions downstream the 
shock, Xi may become much smaller than 1. 

As noted, the magnetic field should have some level of 
inhomogeneity: the scattering centers that isotropise parti- 
cles in both sides of the shock. Magnetic inhomogeneities 
can develop upstream due to relativistic particle streaming 
(e.g. Lucek & Bell 2000), but they should not be suppressed 
by wave damping (see Reville et al. 2007). These inhomo- 
geneities would be advected downstream the shock, thus 
they would be also present there. 

For very high densities, Coulombian/ionization energy 
losses should not suppress acceleration from suprathermal 
energies (e.g. Drury et al. 1996). For that, the following- 
condition should be provided: 

t ion - 3 x 10 11 (mp, e /m e ) -E GcV nZ^ 3 FZq 1 s > i gain , (24) 
which implies: 

n- hC < 4 x 10 9 (m Pje /m e ) FZq 1 v 2 s B_ 3 cm' 3 . (25) 
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The detection of radio emission from electrons with en- 
ergy: 



E 



O.Gvl^Bll^GeV »0.511MeV, 



(26) 



where v is the radiation frequency (1/5 GHz = GHz)), is 
evidence for efficient particle acceleration and hence some 
degree of ionization and B-inhomogeneity at least in some 
sources. However, d, which relates to the i3-inhomogeneity, 
may be small or even energy dependent, not allowing parti- 
cle acceleration to be efficient beyond energies at which elec- 
trons produce synchrotron radio emission. Also, if Z coo i <C 
i?j or mixing were very efficient (large F) , electron acceler- 
ation could also stop at radio emitting energies because 
of large densities and strong relativistic Bremsstrahlung 
losses. In such a case, protons could not reach very high 
energies neither due to strong pp cooling. 

5. Non-thermal emission in the lobes 

5.1. The fate of accelerated particles 

The electrons and protons accelerated by DSA have at in- 
jection an energy spectrum Q(E) cx E~ T up to E maK , with 
r ~ 2 (e.g. Drury 1983) and total luminosity L^f < L s . 
They mainly accumulate downstream the shock, in a re- 
gion that here will be considered homogeneous and with a 
typical size ~ i?j . When protons are present, secondary elec- 
trons and positrons (e^ 1 ) are injected with almost a power- 
law in energy in the region in which these protons interact 
significantly. The injection luminosity and the maximum 
effective energy of e^- will be about a half of the luminos- 
ity going to 7r°-decay gamma-rays and ~ 0.1 the energy of 
protons, respectively (see Kelner et al. 2006). All these par- 
ticles evolve under the downstream magnetic, matter and 
radiation fields, losing energy in the form of synchrotron ra- 
diation, relativistic Bremsstrahlung and IC emission in the 
case of electrons, and high-energy photons, neutrinos, 
and other secondary particles via pp collisions in the case 
of protons. 

Only the radiation from the region closer than i?j 
from the shock will be computed. At distances larger than 
Rj, particles cool through adiabatic losses due to the re- 
expansion of the shocked material, producing scarce emis- 
sion. If densities are large enough (i.e. strongly radia- 
tive shocks, F ^> 1), electrons will cool fast via ioniza- 
tion/Coulombian losses and relativistic Bremsstrahlung, 
and protons through pp collisions, and they will not reach 
the adiabatic cooling region far downstream. 

Far upstream the bow shock, particles with idiff shorter 
than the dominant cooling timescale ti oss could escape into 
the cloud. Some of the highest energy electrons and protons 
would escape from the accelerator in this way. 

The particle distribution, N(E, t), can be obtained solv- 
ing the transport equation (e.g. Ginzburg & Syrovatskii 
1964): 



dN(t,E) d[b(E)N(t,E)} N (t, E) 



dE 



Q{E) 



(27) 



where Q(E) is assumed to be constant in time and b(E) 
includes all the cooling rates —E/t\ oss relevant for N(E,t), 
i.e. synchrotron, relativistic Bremsstrahlung and IC pro- 
cesses for electrons and pp collisions for protons. The adi- 
abatic cooling, which would operate far from the shock 



as described above, is not considered here. The parame- 
ter t esc , the escape time, is the advection timescale (see 
Eq. [6]), different from t^is, which was used to compute 
£ max jiff. This escape time is the dominant timescale for 
the removal of particles from the emitting region. Since in 
general t csc <C inf c , particles will be in the steady state, i.e. 
8N(t,E) /dt = 0. 

5.2. Radiation luminosities and spectral energy distributions 
5.2.1. Luminosities 



Adopting an efficiency 



for transferring 



shock luminosity to non-thermal particles (electrons or pro- 
tons), where f^f < 1, we get: 



f n f L s = 10 35 /„ e £iWrg a" 1 , 



(28) 



where L s3 q — L s /(10 36 ergs -1 ). Since particle escape has 
the dominant dynamical timescale, we can roughly calcu- 
late the relativistic Bremsstrahlung and pp photon bolo- 
metric luminosities in the lobe: 

L Blems ~ 10 33 /„ t0 ! Flo Rjw v & n i,c3 L s36 erg s _1 , (29) 



L pp ~ 10 32 f% w l Fio Rjie v^ 1 n hc3 L s36 erg s 1 



(30) 



These luminosities cannot be higher than L^ t and ~ 1/3 L^ t 
for relativistic Bremsstrahlung and pp collisions, respec- 
tively. These two mechanisms contribute mainly to the 
high-energy part of the spectrum. Regarding the syn- 
chrotron/IC luminosities, in general, when 

f(£max) = tsync,Ic/min[tesc)*BremS)*IC,sync] < 1 : (31) 

Lsyncjc ~ 10 35 /^o.i L s3e erg s" 1 . (32) 
Otherwise: 

(33) 



£ S ync,IC ~ 10 35 /nto.i r 1 L s3 6 erg s 



The IC luminosities will be a minor component unless 
w_9 3> 1. Synchrotron and IC contribute to the low- and 
the high-energy parts of the spectrum, respectively. Note 
that and /„ t may actually be very different in some 
sources (as inferred from the proton-to-electron number ra- 
tio in cosmic rays, i.e. o ~ 100, Ginzburg & Syrovatskii 
1964). 

The highest energy electrons and protons may escape 
from the lobe and radiate in the cloud, although the cor- 
responding luminosities depend on the escape probability, 
which is difficult to quantify. 

5.2.2. Spectral energy distributions 

If synchrotron or IC losses dominate at E < E max for elec- 
trons, there is a break in the particle energy distribution 
N(E,t) at E ~ Eb, in which i S y nc ,ic becomes the shortest 
timescale. Above E^, N(E,t) cx E~( r+1 \ which yields a 
spectral energy distribution (SED) of the radiation that is 
cx e( 2 ~ r )/ 2 for synchrotron and IC dominance, and cx e 1_r 
for relativistic Bremsstrahlung (e is the photon energy). 

For electron energies E < E^, and at any energy for 
protons, advection, and relativistic Bremsstrahlung or pp 
collisions, lead to N(E,t) cx E~ r , which yields a SED 
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oc e 



(3-r)/2 



(synchrotron/IC), and oc e 2 r (both relativis- 6. Thermal emission in the lobes 



tic Bremsstrahlung and pp collisions). Below e 



0.5 MeV and e 



140 MeV, the SEDs of relativistic 



Bremsstrahlung and pp collisions become roughly oc e and 
oc e 2 , respectively. 

Dominant ionization/Coulombian losses, relevant only 
for electrons in our context, lead to N(E, t) oc E 1 ^ , yield- 
ing a SED oc e 2 ~ r / 2 and oc e 3_r for synchrotron/IC and 
relativistic Bremsstrahlung, respectively. 



5.2.3. Requirements from observed spectra 

The non-thermal nature of the observed radio spectra in 
several specific cases implies that strong free-free absorp- 
tion of radio emission should not occur in the lobe or in 
the surroundings. Far from the lobe the ionization degree 
should be small. Close to the lobe, the medium is ionized 
and it is necessary to account for free- free absorption, which 
is expected to be large downstream the bow shock. The free- 
free opacity there can be written as (Rybicki & Lightman 
1979): 



T ff - 0.03 T, 



-1.5 -2 



( Hz ^iO.l n c5 ^10 'l6 ) 



(34) 



where T 5 = T/10 5 K, and l 16 = l/(10 16 cm) is the typical 
size of the region. Also, suppression of emission at frequen- 
cies v < 2 x 10 8 rLj^Xio.i/ B_3 Hz, due to the Tsytovich- 
Razin effect, should be considered. 

Another condition that should be fulfilled given 
the observed radio spectra is that either electron es- 
cape or relativistic Bremsstrahlung should dominate over 
Coulombian/ionization losses at low electron energies, i.e. 
a ~ 0.5 (F v oc v~ a ). From the timescale ratios: 



tu 



iAcsc ~ 8 x 10 2 E GcV n U c 3 -Flo 1 v i-,Ls and 



Brcms 



9E, 



GoV : 



(35) 
(36) 



it is seen that upstream and downstream the shock, ra- 
dio emitting electrons are dominated by escape or relativis- 
tic Bremsstrahlung losses. Note that sources with spectra 
harder than F v oc v 



-0.5 



may be still explained in a non- 
thermal scenario by moderate free- free absorption and/or 
ionization cooling (or a thermal component, see below). 



5.3. Deriving the magnetic field strength 

Assuming a value for L„ t , plus some additional simplifying 
assumptions, a formula has been obtained to derive the 
magnetic field strength consistent with the observed radio 
fluxes. Taking the radio fluxes at a certain frequency, a 
particle energy distribution with T ~ 2 and normalized 
with the total energy ~ L^ t t csc , the synchrotron power for 
one electron (E ps 4.1 x 10~ 15 B 2 _ 3 Eq cV erg s" 1 ), and the 
reasonable simplification that electrons of energy E produce 



only photons of frequency v 
obtain: 



5 x 10 9 B_ 3 El 



GcV 



Hz. 



B 



0.04(^ t35 )- 2 / 3 i? j 



2/3 2/3 1/3 
16 V " 



j4/3 



^5 GHz ^3 kpc 



F^ Jy mG,(37) 



where <ikpc is the distance the source in kpc, and F v m j y 
the flux in mjy at the relevant frequency. In case the ra- 
dio emitting leptons are secondary e ± , then in Eq. (|37|) it 
should be changed: 



0-1 (t esc /tpp) 



(38) 



The shocked material is heated up to T « 2.3 x 
10 7 w 2 g K (2.4 1> 2 8 keV) and generates line and thermal 
Bremsstrahlung continuum emission. The total luminosity 
cannot overcome L Sl and it might be much smaller if the re- 
verse shock is adiabatic and the bow shock relatively slow. 
This radiation can be absorbed in the cloud core, since the 
photo-electric opacity coefficient is rx ~ 20(ek e v)~ 2 ' 5 -^H23 
(within a factor of 2 in the relevant energy range; see 
Morrison & McCammon 1983), where ekcV = e/keV is the 
photon energy and iV H = 10 23 iV H 23 = 10 23 n c5 l 18 cm -2 the 
neutral hydrogen column density of the lobe surroundings. 

The density in the shock regions should neither be too 
high, to avoid suppression of the acceleration due to strong 
cooling, nor too low, to avoid that va iS V Tt bs- Actually, 
there is room for the shocks to be radiative, producing ther- 
mal X-rays that may escape the cloud, without necessarily 
suppressing particle acceleration. 

The fact that thermal radio emission should not over- 
come the non-thermal component at the same frequencies 
put some constraints on the scenario. The thermal SED 
peaks at energies up to ~ keV and is harder than the 
synchrotron SED, which should generally peak at lower 
energies. This means that the total L S ync must be <C L s 
(-L S ync ^ 10 -4 L s for r = 2), if the non-thermal radio emis- 
sion is to be dominated by the thermal one. This condition 
is hard to fulfill unless f£ t <C 1. It is worth noting that a 
radio spectra harder than z/ -0 - 5 may be non-thermal radi- 
ation contaminated by a thermal component. 

As mentioned in Sect. [TJ several massive YSOs present 
non-thermal radio emission. Two among the most relevant 
of them are studied in next section: IRAS 16547—4247 and 
the complex source HH 80—81. 



7. Application to IRAS 16547-4247 North and 
HH 80 

The model described in Sect. 15.11 is applied now to the 
northern radio lobe of the massive YSO IRAS 16547-4247 
(IRAS-N), and to the radio lobe HH 80 in the complex 
source HH 80—81. Both lobes have a clear non-thermal 
nature (e.g. GAR03; MRR93). IRAS-N shows also an ex- 
tended structure pointing to the South-East, and HH 80 has 
a similar non-thermal source very nearby, HH 81. We will 
not consider here either the South-East extension of IRAS- 
N nor HH 81. Note however that the extension in IRAS-N 
may be in fact a fore/background object (RMF08), and 
HH 81 could be part of the expected complex structure of 
the jet termination region (see Sect [3j see also Heathcote 
et al. 1998 -HRR98- for a detailed optical study of the 
HH 80-81 complex). 

7.1. IRAS-N and HH 80 properties 

We focus on IRAS-N and HH 80 because both sources are at 
the two extremes of the density parameter range presented 
above. IRAS-N is embedded in a very dense cloud, with 
n c ~5x 10 5 cm" 3 (GAR03), whereas HH 80 is thought to 
be close to the border of a cloud, in a more diluted medium 
with n c - 10 2 - 10 3 cm" 3 (e.g. MRR93; HRR98; Pravdo 
et al. 2004 -PTM04-). The distances to IRAS-N and HH 80 
are d - 2.9 kpc and - 1.7 kpc (RMF08; MRR93), and 
the central stars show luminosities of w 5 x 10 38 and 
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8 x 10 37 erg s" 1 , respectively (GAR03; MRR93). This radi- 
ation provides the main contribution to the infrared emis- 
sion in the lobes, yielding photon energy densities there of 
u ~ 2 x 10" 9 erg cm" 3 for IRAS-N, and 2 x 10" 12 erg cm" 3 
for HH 80. The distances from the central star to the lobes 
are Zj w 5 x 10 17 cm for IRAS-N, accounting for a jet incli- 
nation angle of 84° (Garay et al. 2007), and ~ 10 19 cm for 
HH 80, with a not so well constrained jet inclination (see 
however HRR98). The lobe sizes for IRAS-N and HH 80 
are about Rj » 1.6 x 10 16 cm and 5 x 10 16 cm, (RMF08; 
MRR93), the velocities of the jets are expected to be around 
Vj ~ 5 x 10 7 cm s" 1 and ~ 10 8 cm s" 1 , and the ve- 
locities of the bow shocks would be Vbs < 10 7 cm s" 1 
and - 5 x 10 7 cm s"\ respectively (RMF08; MRR93, 
MRR95, HRR98). IRAS-N has not been detected in X-rays 
(see ARA07), whereas HH 80 has been detected by XMM 
(PTM04). 

7.1.1. Derived parameters 

The values of v hs for IRAS-N imply that \ & 0.04 and 
u r ~ v- y In fact, from the inferred age of IRAS-N, tm c ~ 
10 11 s (Garay et al. 2007), Zj and Vj, a particle density 
of rij ~ 5 x 10 3 cm" 3 can be derived, i.e. \ ~ 0-01 and 
«bs ~ 5 x 10 6 cm s" 1 , consistent with the limit given above. 
Such a value for rij, together with Rj and Vj, renders a 
Lj ~ 5 x 10 35 erg s" 1 for IRAS-N. Given that v x ~ uj in 
this source, it is the case that L s ~ Lj in the reverse shock. 
The bow-shock luminosity will be ~ 4 x 10 34 erg s . 

For HH 80, since n c ~ nj, then % ~ 1, and therefore we 
have that v x ~ v\, s . In the case of HH 80 this means that 
both the reverse shock and the bow shock may contribute 
to the non-thermal radiation (see also HRR98). Accounting 
for x an d fj, it can be inferred ixife ~ 3 x 10 11 s, not far 
from the value discussed in MRR93; also, v r ~ v\, s . Taking 
Lj ~ 2 x 10 36 erg s" 1 (MRR95), we get nj - 4 x 10 2 cm" 3 
and, since x ~ 1, n c ~ 4 x 10 2 cm" 3 . This latter value is 
between those inferred from X-ray and optical observations 
(see the discussion in PTM04). We will treat the bow shock 
and the reverse shock in HH 80 as a single physical system, 
with speed 5 x 10 7 cm s" 1 and L s — L y The full list of the 
relevant properties of IRAS-N and HH 80, together with 
the derived parameters, is presented in Table [TJ 

Interestingly the central star is brighter in IRAS-N than 
in HH 80, but Lj seems smaller in the former. This could be 
related to the larger density of the environment in IRAS- 
N. This may have induced jet deceleration through, e.g., 
medium mass entrainment in the jet. Nevertheless, the un- 
certainties are large and no strong conclusions can be de- 
rived in this regard. 

It is worth mentioning that very powerful, slow molec- 
ular outflows with luminosities 10-100 times larger than 
I/j have been detected in the two sources (see RMF08 and 
references therein). 

7.2. Estimates of the emission in IRAS-N and HH 80 

7.2.1. Constraints on the non-thermal population 

In order to model the radio emission from IRAS-N and 
HH 80, and to compute the radiation at high and very 
high energies, the values of B and L^f f° r both sources 
are required. B should be well below B eq (see Eq. [5J. The 



Table 1. IRAS-N and HH 80 properties and derived pa- 
rameters (see the text for details) 







xlxl oU 


n c [cm 3 ] 


5 x 10 5 


4 x 10 2 


d [kpc] 


2.9 


1.7 


L* [erg s" 1 ] 


5 x 10 38 


8 x 10 37 


u [erg cm -3 ] 


2 x 10~ 9 


2 x 10" 12 


Zj [cm] 


5 x 10 17 


1Q 19 


U l l C1T1 J 


C v 1(1 5 


a v i n 2 


Rj [cm] 


1.6 x 10 16 


5 x 10 16 


Vj [cm s" 1 ] 


5 x 10 7 


10 s 


v ba [cm s" 1 ] 


5 x 10 6 


5 x 10 7 


Vi [cm s" 1 ] 


5 x 10 7 


5 x 10 7 


tlife [s] 


10 11 


3 x 10 11 


Uj [cm" 3 ] 


5 x 10 3 


4 x 10 2 


X 


0.01 


1 


Lj [erg s" 1 ] 


5 x 10 35 


2 x 10 36 



non-thermal luminosity L^jf will be taken equal to 0.1 L s , 
or i = l- The observed radio fluxes and spectral indices 
are ~ 8.7 and 3 mJy at 8 GHz and a ~ 0.5 and 0.3, respec- 
tively (GAR03; MRR93). The value of the index T for the 
radio emitting particles can be obtained from a, as shown 
in Sect. 15.11 Fixing L^t an d knowing the radio fluxes, B 
can be estimated to zeroth order with Eq. (|3"T)) . We note 
that equipartition with the relativistic particles would lead 
to magnetic fields of the order of 0.1-1 mG in the emitting 
regions (see ARA07 for the case of IRAS-N). 

The nature of the radio emitting particles may be pri- 
mary electrons accelerated in the shock, or secondary 
from pp collisions. If primary electrons dominate the pro- 
duction of the radio emission, the relativistic proton popu- 
lation will be constrained by the fact that secondary radi- 
ation cannot overcome that of primary electrons (although 
pp collisions may still be a significant source of high en- 
ergy emission) . If secondary were the origin of the radio 
emission, the value of V for protons should be slightly softer 
than that of secondary e^ 1 (Kelner et al. 2006), and the ra- 
dio contribution from primary electrons should be minor. 

7.2.2. The emission in IRAS-N 

Given the values of n c and i>b s in IRAS-N, it seems un- 
likely that the bow shock is accelerating these electrons. 
This shock will be strongly radiative, peaking the ther- 
mal emission in the optical/UV, with a luminosity of few 
times 10 34 erg s" 1 . On the other hand, the reverse shock is 
marginally radiative because / coo i ~ Rj, peaking at 0.5 keV 
with a bolometric luminosity ~ 10 35 erg s" 1 . As men- 
tioned above, the photoelectric absorption rx is very large, 
likely > 100, implying an almost complete suppression of X- 
rays, which would explain the non-detection of the source. 
Concerning the non-thermal radiation, from /^J 1 = 1 one 
obtains L e at p - 5 x 10 34 erg s" 1 . 

If primary electrons produced the radio emission, the 
magnetic field should be B ~ 0.1 mG for the adopted 
/ n t, and the maximum energy, limited by diffusive escape, 
B max ~ 4 x 10 2 GeV for both electrons and protons. If 
secondary produced the radio emission, B > 0.5 mG, 
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and -©max > 2 x 10 3 GeV (only protons), limited again by 
diffusive escape. 

The density downstream the bow shock may be up to 
10 4 times larger than downstream the reverse shock due 
to strong radiative cooling in the former region. Since £rt 
should be much shorter than tnt e , we assume that the RT 
instabilities have time to develop mixing reverse and bow 
shock downstream material, effectively increasing the den- 
sity in the reverse shock. We have adopted a value for F 
such as that the resulting luminosities due to relativistic 
Bremsstrahlung and pp collisions be of the order of L nt . In 
this way, F is used as a free parameter that is optimized to 
get high gamma-ray luminosities, not to suppress accelera- 
tion, and accounting for the constraints given by the obser- 
vational data at lower energies and the F-limits provided 
in Sect. 13.21 Following such an approach we have adopted 
F ~ 25. Note that, despite the high density, the Tsytovich- 
Razin effect may be neglected, since the entrained dense 
bow shock material should be cold and hardly fully ion- 
ized. 



7.2.3. The emission in HH 80 

In the case of HH 80, both the reverse and the bow shock are 
adiabatic. We predict thermal X-rays peaking at ~ 0.5 keV 
with intrinsic luminosities of the order ~ (i?j/Z coo i) L s ~ 
10 34 erg s _1 . This value is ~ 100 times bigger than that 
given by PTM04, where no intrinsic absorption was as- 
sumed. However, if an intrinsic jVh ~ 5 x 10 21 cm~ 2 in 
the surroundings of HH 80 were adopted, the factor of 
100 could be explained with photoelectric absorption (see 
the strong dependence of ^ coo i and Tx on photon energy 
and w r .bs)- To compute the non-thermal emission, as noted, 
the two shocks are treated as just one with ftg x = 1, i.e 
C~ 2 x 10 35 ergs" 1 . 

In the case dominated by primary particles, the mag- 
netic field would be B ~ 0.005 mG for the adopted 
/ n t, with a maximum energy limited by diffusive escape 
E mSLX ~ 60 GeV for both electrons and protons. In the 
case of dominance by secondaries, B > 0.02 mG, and 
-Emax > 3 x 10 2 GeV (only protons), limited again by diffu- 
sive escape. Since the medium is quite diluted, relativistic 
Bremsstrahlung and pp collisions are not as efficient as in 
IRAS-N, but the fact that and t esc are both larger 
renders not so different values for Lbrems.pp- 



7.3. Computed spectral energy distributions 

In Figs. H EE OS and [SJ the SEDs computed for IRAS-N 
and HH 80 are shown. Two scenarios are adopted for both 
sources, one in which the radio emission is dominated by 
primary electrons, and another one in which the dominant 
radio emitters are secondary e ± . In the former, L^ t and 
L^ t has been taken equal yielding a ~ 1; in the latter, 
we have derived just a lower-limit for a to avoid primary 
emission to be significant. In IRAS-N, if secondary e ± were 
the source of radio emission, a should be > 10. In HH 80, 
given the relatively low densities and high magnetic fields 
of the secondary e ± scenario, the proton to electron number 
ratio a should be > 1000. 

It is remarkable that the high-energy components in the 
SEDs, associated with relativistic Bremsstrahlung and/or 
pp collisions, have significant luminosities in the high- 
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Fig. 2. Spectral energy distribution of the non-thermal 
emission for IRAS-N in the primary electron scenario. 
The IC contribution is negligible and not shown here. 
Observational points are from IRAS 16547—4247 (radio, 
Rodriguez et al. 2005; X-rays, ARA07). The 1 yr/5 a sen- 
sitivity of Fermi in the direction of the galactic plane is 
shown. A curve above 100 GeV showing a luminosity cor- 
responding to 0.01 Crab, typical sensitivity of a Cerenkov 
telescope for exposures of ~ 50 hr, is also presented. 
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Fig. 3. The same as in Fig.[5]but for the secondary e ± case. 



energy and very high-energy range and fulfill the X-ray 
constraints. The synchrotron emission peaks in the opti- 
cal/UV, and can be the dominant cooling channel of elec- 
trons only if the magnetic field is rather high and densi- 
ties low. Interestingly, in the primary electron scenario of 
HH 80, the electron component does not achieve energies 
beyond those to emit radio synchrotron emission, but rel- 
ativistic Bremsstrahlung and pp emission may still be sig- 
nificant at GeV energies. 

The list of the parameter values adopted to calculate the 
SEDs, together with the radio properties of the sources, 
is presented in Table 1. The parameter values have been 
adjusted numerically and are slightly different from those 
given in Sect. 17.21 



8. Detectability 

The SEDs presented in Figs. 0(5] show that massive YSOs 
can produce significant amounts of gamma rays, although 
the results are quite sensitive to the available densities 
downstream. For sources with high densities like IRAS-N, 
if the development of RT instabilities does not interfere 
with particle acceleration, relativistic Bremsstrahlung and 
pp collisions will be quite efficient. For low-density sources 
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IRAS-N 


HH 80 


a 

F s GHz [mjy] 


0.5 
8.7 


0.3 
3 




primaries" secondaries 6 


primaries secondaries 



a 


l 


> 10 


1 


> 1000 


n [cm -3 ] 


5 x 10 5 


5 x 10 5 


1.6 x 10 3 


1.6 x 10 3 


^esc [s] 


9 x 10 8 


9 x 10 8 


4 x 10 9 


4 x 10 11 


L e nt p [erg s- 1 ] 


5 x 10 34 


5 x 10 34 


2 x 10 35 


2 x 10 35 


B [mG] 


0.25 


2 


0.003 


2.5 


£ max (e) [GeV] 


7 x 10 2 


~ 5 x 10 2 


35 


~ 10 3 


£max (p) [GeV] 


7 x 10 2 


5 x 10 3 


35 


10 4 


r 


2.2 


2.4 


1.6 


1.8 



a Primary electron scenario. 
b Secondary e scenario. 
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Fig. 4. The same as in Fig.f2]but for HH 80. Observational 
points are from Marti et al. (1993) (radio) and PTM04 (X- 
rays). The X-ray detected point is shown as an upper- limit. 
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Fig. 5. The same as in Fig. [3] but for HH 80. 



accordance with estimates derived through Zeeman mea- 
surements (e.g. Crutcher 1999). In the case of low-density, 
HH 80- like, sources, the situation may be more complicated. 

The primary electron scenario of HH 80 requires a very 
low magnetic field in the shocked regions. As noted, low 
densities imply large lobes and therefore lower magnetic 
fields in the jet head. However, the value of B in the bow 
shock can hardly be smaller than that in the cloud, which 
is expected to be, given the cloud densities, several times 
higher than the one adopted in our calculations. Therefore, 
in such a scenario and source type, if detected, gamma rays 
would likely come from the reverse shock. Otherwise, in the 
secondary e ± scenario in HH 80 the magnetic field must be 
quite high, ~ 2.5 mG, regardless the shock involved, below 
but close to the maximum value (see Eq. [8]). Furthermore, 
the value for a required in this case, ii 1000, may be too 
large. Also, the hard particle energy distribution required 
may be difficult to explain in the context of linear theory of 
Fermi I particle acceleration. Such a hard radio spectrum 
may be explained by marginal free-free absorption or by 
an additional thermal component, but then the expected 
non-thermal fluxes at higher energies would be smaller due 
to a softer particle energy distribution. In any case, despite 
these caveats, one cannot still rule out HH 80 and similar 
objects as gamma-ray emitters. 

It is worth noting that the assumptions adopted in this 
work are quite conservative. The parameter uncertainties 
are relatively large, and a more optimistic, but yet consis- 
tent with observations, choice of densities, shock velocities 
and jet luminosities, could easily move the SED curves up 
by a factor of several. 



like HH 80, the efficiency of relativistic Bremsstrahlung and 
pp collisions is lower, but the lobe sizes are expected to be 
large, increasing the escape timescales and thereby the ra- 
diation outcome (see Sect. I5.2p . Therefore, for f^ p 0.1, 
it can be expected that the termination regions of massive 
YSO jets will be eventually detected by Fermi and also by 
Cherenkov telescopes through long enough exposures. 

In our calculations, the magnetic field strength B and 
Ljjjf have been adjusted for both to explain the radio fluxes 
and to obtain significant gamma-ray fluxes. The magnetic 
field assumed for high-density, IRAS-N-like sources is in 



9. Discussion 

Romero (2008) pointed out that the detection of massive 
protostars at gamma-ray energies would open a new win- 
dow to star formation studies. The detection of the cutoff 
in the SED would give important insights on the accel- 
eration efficiency in the terminal shocks of the outflows. 
The SED can also shed light on the densities, magnetic 
fields, velocities, and diffusion coefficients in the shocked 
regions. Although we do not expect that massive proto- 
stars should be among the bright sources detected by Fermi 
(Abdo et al. 2009), our calculations show that they could 
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show up in further analysis of weaker sources after few 
years of observation. The emission levels above 100 GeV, 
close to 0.01 Crab, could be detectable by current and fu- 
ture Cherenkov telescopes for observation times moderately 
longer than 50 hours. 

However, not only massive protostars, but also the re- 
gions in which they form, may be gamma-ray emitters. As 
mentioned in Sect. [H some amount of the highest energy 
particles may escape to the cloud far upstream the bow 
shock. It is hard to estimate the fraction of electrons and 
protons that would be released in the cloud, which depends 
strongly on the diffusion coefficient of the pre-shock cloud 
medium and the bow shock velocity and size. However, they 
might carry a non-negligible fraction of L^f if T ~ 2. In 
fact, in the case that L„jf ~ 0.01 Lj were in the cloud, 
massive YSOs may inject an amount of protons well above 
the average galactic level at several hundreds of GeV, and 
the radiation resulting from pp may be detectable (for a 
general case, see Aharonian & Atoyan 1996), competing 
with that produced in the lobe itself. For leptons, the emis- 
sion at high energies may be relevant for low magnetic 
fields, i.e. when the maximum energy is determined by dif- 
fusive escape and dominant relativistic Bremsstrahlung in 
the cloud. The spectrum of the gamma rays, generated by 
pp collisions for protons and relativistic Bremsstrahlung for 
electrons, should be very hard since only the highest en- 
ergy particles escape, peaking at e ^ E ma ^. The cloud syn- 
chrotron emission should be quite diluted and dominated 
by the lobe. 

A clustering of gamma-ray sources should be present 
in regions with large molecular clouds and star formation, 
as already inferred from EGRET data (e.g. Romero et al. 
1999). The accumulation of cosmic rays accelerated in the 
radio lobes into the molecular cloud can produce extended 
gamma-ray sources. These radio lobes may be difficult to 
detect. Neither UV nor hard X-ray counterparts related to 
thermal Bremsstrahlung produced in the shock downstream 
regions are expected to be observed from these sources be- 
cause of the large absorption and/or low emission levels. 
Deep inside the cloud, even radio emission may be missing 
due to strong free-free absorption, so the exact number of 
accelerators could be hard to estimate. Also, cosmic-ray re- 
acceleration inside the clouds due to magnetic turbulence 
(e.g. Dogiel et al. 2004) could result in stronger sources. 
Therefore, the combined effect of several protostars deeply 
embedded in giant clouds might be responsible for GeV- 
TeV sources found in star forming regions by EGRET, 
Fermi, AGILE and Cherenkov telescopes. We conclude that 
massive clouds with high IR luminosities and maser emis- 
sion (tracers of massive star formation) deserve detailed 
study with Fermi and ground-based Cherenkov telescopes. 
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